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Abstract 

(N 

^ ' We present a general formalism that provides a systematic com- 

^^ ■ putation of the linear and non-linear perturbations for an arbitrary 

T-|- ! number of cosmological fluids in the early Universe going through 

VO I various transitions, in particular the decay of some species (such as 

t"~~. ■ a curvaton or a modulus). Using this formalism, we revisit the ques- 

^P , tion of isocurvature non-Gaussianities in the mixed inflaton-curvaton 

scenario and show that one can obtain significant non-Gaussianities 
dominated by the isocurvature mode while satisfying the present con- 
straints on the isocurvature contribution in the observed power spec- 
trum. We also study two-curvaton scenarios, taking into account the 
C^ ' production of dark matter, and investigate in which cases significant 

non-Gaussianities can be produced. 

1 Introduction 

In many occasions, cosmology has been and still is an invaluable means 
to constrain particle physics models. These constraints can arise by using 
information from homogeneous cosmology, such as the constraints on the 
light degrees of freedom at nucleosynthesis. With the discovery of the Cosmic 
Microwave Background (CMB) fluctuations, new constraints arise from the 
observed power spectrum of linear perturbations. Even more recently, the 



O 



>$ 



upper bounds on primordial non-Gaussianities have started to be used to 
constrain early Universe scenarios. 

Although the simplest early Universe models are based on inflationary 
models with a single scalar field, many models involve additional scalar fields, 
which can play a dynamical role during infiation or simply be spectactor fields 
(see e.g. [1] for introductory lectures). The existence of several degrees of 
freedom opens up the possibility of isocurvature perturbations, i.e. perturba- 
tions in the particle density ratio between two fiuids, for example cold dark 
matter (CDM) isocurvature perturbations (between CDM and radiation) or 
baryon isocurvature perturbations (between baryons and radiation). Since 
primordial isocurvature perturbations leave distinctive features in the CMB 
anisotropies, they can be in principle disentangled from the usual adiabatic 
mode. The present upper bound on the isocurvature contribution to the 
power spectrum provides a stringent constraint. 

This is the case for the curvaton scenario [2] where large residual isocur- 
vature perturbations (for CDM or baryons) can be generated, depending on 
how and when CDM or baryons are produced [3, 4] (see also [5, 6] for more 
detailed scenarios). The same constraints apply to moduli that are light 
during infiation, and thus acquire super-Hubble fiuctuations, as discussed 
recently in [7]. 

Another potentially useful information on primordial perturbations is the 
amplitude and shape of their non-Gaussianity. So far, the current CMB data 
seem to favour a non-zero amount of so-called local non-Gaussianity [8] , but 
Planck data will be needed to confirm or infirm this trend. Several models 
can generate local non-Gaussianity (see e.g. [9] for a recent review): multiple 
field infiation (during infiation or at the end of infiation: see e.g. [10]), 
modulated reheating [11, 12], curvaton, modulated trapping [13], etc. It is 
thus interesting to combine the constraints on isocurvature modes and non- 
Gaussianity to explore the early Universe physics, as has been done recently 
in various scenarios [14, 15, 16, 17, 18, 19, 20, 21]. 

The purpose of the present work is to give a unified treatment of linear and 
nonlinear perturbations, which enables to compute their evolution through 
one or several cosmological transitions, such as the decay of some particle 
species. Our treatment takes into account the various decay products and 
their branching ratio. Our formalism can thus be applied to a large class of 
early Universe scenarios, in order to compute automatically their predictions 
for adiabatic and isocurvature perturbations, and their non-Gaussianities. 
As input, one simply needs parameters that depend on the homogeneous 
evolution. This thus provides a simple way to confront an early Universe 
scenario, and its underlying particle physics model, with the present and 
future cosmological data. 



As applications to our general formalism, we consider two specific ex- 
amples. The first example is a more refined treatment of the isocurvature 
perturbations and their non-Gaussianity in the mixed curvaton-infiation sce- 
nario [22, 23, 24]. The second example deals with a multiple-curvaton sce- 
nario [25, 26, 27, 28]. In both examples, we generalize the results that have 
been obtained in previous works, allowing the curvaton to decay into several 
species. 

This paper is organized as follows. In Section 2, we introduce the non- 
linear curvature and isocurvature perturbations. Section 3 is devoted to the 
general treatment of a cosmological transition, such as the decay of some 
particle species. In Section 4, we focus on the first application, namely the 
mixed curvaton-infiaton scenario with a single curvaton. In Section 5, we 
consider scenarios with two curvatons. We conclude in the final Section. 

2 Non-linear curvature perturbations 

We first introduce the notion of non-linear curvature perturbation. Several 
definitions have been proposed, which turn out to be equivalent on large 
scales, and we will follow here the covariant approach introduced in [29, 30], 
and reviewed recently in [31]. 

For a perfect fluid characterized by the energy density p, the pressure 
P and the four- velocity -u", the conservation law V oT\ = for the energy- 
momentum tensor. Tab = (p + -P) UaUb + Pgab, implies that the covector 

N 

C^ = V^N--VaP (1) 

P 

satisfies the relation 

with the definitions 

e = V,M^ N = ^fdrQ, (3) 

where r is the proper time along the fluid worldlines and a dot denotes a Lie 
derivative along -u", which is equivalent to an ordinary derivative for scalar 
quantities (e.g. p = u°''VaP)- N can be interpreted as the number of e-folds 
of the local scale factor associated with an observer following the fluid. 

The covector (a can be defined for the global cosmological fluid or for 
any of the individual cosmological fluids, as long as they are non-interacting 



(the case of interacting fluids is discussed in [32]). Using the non-hnear 
conservation equation 

P = -3iV(p + P) , (4) 

which follows from u^VaTI = 0, one can re-express (a in the form 

li w = P/p is constant, the above covector is a total gradient and can be 
written as 

C = V, N+ /\ Inp . (6) 

On scales larger than the Hubble radius, the above definitions are equiv- 
alent to the non-linear curvature perturbation on uniform density hypersur- 
faces as defined in [33], 

J-p P ^Jp (1+W)p 

where H = a/a is the Hubble parameter. 

It will be useful to distinguish the non-linear curvature perturbation ( of 
the total fiuid, from the individual non-linear perturbation (a that describes 
the cosmological fiuid A (with wa = Pa/ Pa = for a pressureless fiuid or 
wa = 1/3 for a relativistic fiuid), defined by 

1 T.. ( Pa 



r=5N+ — In , 

3(l + ti;J Vp. 

where a bar denotes a homogeneous quantity. 

Inverting this relation yields the expression of the inhomogeneous energy 
density as a function of the background energy density and of the curvature 
perturbation ^^, 

P.=P.e^(^+--)«--^^), (9) 

which we will use many times in the following. 

The non-linear isocurvature (or entropy) perturbation between two fiuids 
A and B is defined by 

Sa,b = 3{U-Cb). (10) 

In the following, we will always define the isocurvature perturbations with 
respect to the radiation fiuid, so that our definition for the isocurvature 
perturbation of the fiuid A will be 

5^ = 3(a-Cr), (11) 

where C,^ is the uniform- density curvature perturbation of the radiation fiuid. 



3 Decay 

Let us now consider a cosmological transition associated with the decay of 
some species of particles, which we will call a. 

In the sudden decay approximation, the decay takes place on the hyper- 
surface characterized by the condition 

H^ = T„. (12) 

Therefore, since H depends only on the total energy density, the decay hy- 
persurface is a hypersurface of uniform total energy density, with SN^^ = (, 
where ( is the global curvature perturbation. Using (9), the equality between 
the total energy densities, respectively before and after the decay, thus reads 

J2 P^-e^(^+--)(^--"'^) = Pdccay = Yl PB+e^^i+--)(^-+-«, (13) 

A B 

where the subscripts — and + label quantities defined, respectively, before 
and after the transition. 

3.1 Before the decay 

The first equality in (13) leads to 

^fiAe='(i+-^)(^^-« = l, (14) 

A 

where we have defined Qa = Pa~ /Pdcca.y (to avoid confusion, the Qa are 
always defined just before the decay). The above relation determines (^ as a 
function of the (a-- 

At linear order, this gives 

C = i$^^ACA- (first order) (15) 

with the notation 

nA = {i + WA)^A, n = J2^A. (16) 



Expanding (14) up to second order, one finds 



c4E"- 



'A 



Ca- + 1{1+wa){Ca--C? 



(second order) (17) 



where, on the right hand side, ( is to be replaced by its first order expression 
(15). 



3.2 After the decay 

We now consider the outcome of the decay. In general, the species a de- 
cays into various species A, with respective decay widths Tacj- Defining the 
relative branching ratios 

7Aa = -f^, r^ = ^ r^^ , (18) 

one can write the energy density of the fiuid A after the decay in terms of 
the energy densities of A and of a as 

PA+ = Pa- + lAuPu- ■ (19) 

Using (9), one can rewrite this nonlinear equation in terms of the curvature 
perturbations Ca+, Ca~ and Co--, which yields 

n. p3(1+«;a)(Ca--C) _L^. 7. p3(l+«><,)(C<.--C) 
g3(l+«;A)(CA+-C) = P'^^^ +lAapa~e ^ 

Pa- + lAaPa- 

This expression thus gives Ca+ as a function of Ca--, Co- and of the global C- 
Substituting C, in terms of Co and of all the Cb- , one finally obtains Ca+ as a 
function of all the Cb-- 

Following this procedure, one finds that the linear curvature perturbation 
for any given fiuid A is given by 

C^+ = Y. ^a'^Cb- (first order) (21) 



with 



B 



A T f , f i'^A - Wa)^A 



^ ii + WA)n ^ ' 

^' = ^'Y^^A^^' {i + WA)n ^''^ 

T/ = / /"-^'"-^^^ C^A- (24) 

In the above expressions, we have introduced the parameter 

which represents the fraction of the fiuid A that has been created by the 
decay. If A does not belong to the decay products of a, then Ja = 0. The 



opposite limit, /^ = 1, occurs when all the fluid A is produced by the decay. 
For the intermediate values of J'a, part of A is produced by the decay while 
the other part is preexistent. 

In the following, we will assume that the decaying species behaves like 
non-relativistic matter (this is the case for a curvaton or modulus fleld that 
oscillates in a quadratic potential) and we will thus always use Wa = 0- 

From the above expressions (22-24), it is straightforward to check that 



E^. 



B 



1. 



(26) 



B 



The post-decay perturbation (a+ can thus be seen as the barycenter of the 
pre-decay perturbations (b- with the weights T^^ (all these coefficients sat- 
isfy < T^^ < 1 for Wo- = 0). Note that if the fluid A is not produced in the 
decay (i.e. /a = 0), then the transfer coefficients are trivial: T^^ = 6j^. 

Since it is convenient to use the same range of species indices before and 
after the transition, we also introduce the coefficients T^^ = 0, which imply 
that Ccr+ = 0. This convention will be especially useful when one needs to 
combine several transitions, as we will discuss soon. 

At second order, expanding (20) and substituting the first order expres- 
sion (15) for (, one obtains 







E^/Cb- 



Eu^^Cb-Cc- 



(second order) 



(27) 



B,C 



with 



ttBC 



Q. 



c , 



T4b(1 + wb)Sbc + 2-^(u'A - wb)Tab - (1 + wa)TabTac 



QrQ 



B^l-C 



Q^ 



1+wa-^ Tad(1 + wd) 



D 



(28) 



The change of the various isocurvature perturbations, defined in (11), can 
also be determined by using the above expressions. In particular, at linear 
order, one finds, using the property (26), the simple expression 



^A+ = Yl i^A^" - Tr "") Sb- (first order) . 



(29) 



3.3 Several transitions 

If the early Universe scenario involves several cosmological transitions, for ex- 
ample several particle decays, one can use the above expressions successively 



to determine the final "primordial" perturbations, i.e. the initial conditions 
at the onset of the standard cosmological era. 

For linear perturbations, the expression of the final perturbations as a 
function of the initial ones, is simply given by 

cy^=E^Ai\ T=nT[.] (30) 

B k 

where T is the matricial product of all transfer matrices T[k] , which describe 
the successive transitions. 

The cosmological transitions can result from the decay of some particle 
species but they can be of other types. For example, if CDM consists of 
WIMPs (Weakly Interacting Massive Particles), the freeze-out can be treated 
as a cosmological transition. If radiation is the dominant species at freeze- 
out, then Cc+ = Cr- But, if other species are significant in the energy budget 
of the universe at the time of freeze-out, any entropy perturbation between 
the extra species and radiation will modify the above relation. The presence 
of a pressureless component, like a curvaton, leads to [4] 

at linear order, while the other C^a remain unchanged. The symbol "o"" de- 
notes here the conglomerate of all pressureless matter at the time of freeze- 
out, except of course the CDM species that is freezing out. 



4 Scenario with a single curvaton 

Let us now apply our formalism to a simple scenario with only three initial 
species: radiation (r), CDM (c) and a curvaton (o\ considered in e.g. [34]. 
After the decay of the curvaton, the radiation and CDM perturbations remain 
unchanged and provide the initial conditions for the perturbations at the 
onset of the standard cosmological phase (let us say around T ~ 1 MeV). 

4.1 Perturbations after the decay 

4.1.1 Linear order 

According to the expressions (22-24), the linear transfer matrix T45 is given 
in this case by 



nn /y» /y» /y» 

1 - /c /c 
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T=\ l-/c /c I , X^ = -J^, Xc = -^cXr (32) 



where the order of the species is (r, c, a). This means that the hnear curvature 
perturbations for radiation and for CDM, after the curvaton decay, are given 
respectively by 

Cr+ = (1 ^ ^rj Cr- + ^c Cc- + [^r "~ ^cj Co— I'J'j) 

and 

Cc+ = (l-/c)Cc- + /cC.- (34) 

The entropy perturbation after the decay is thus 

-^c+ = Cc+ ~ Cr+ = [i- ~ Jc~ Xc)Cc~ + y^r ~ IjCr- + [jc + ^^c — Xr)Ca- , (35) 

which can also be expressed directly in terms of the pre-decay entropy per- 
turbations, following (29), 

Sc+ = (1 - /c - Xc) Sc- + ifc + Xc- Xr) S„- . (36) 

Note that, if many CDM particles are created by the decay of the curvaton, 
a significant fraction of them could annihilate, leading to an effective sup- 
pression of the final isocurvature perturbation. This effect has been studied 
in [5] and can easily be incorporated in our formalism. 

In practice, we will need the above expressions only in the limit Xc = 
since fie is usually negligible when the decay occurs. The coefficient Xr, which 
we will shorten into r from now on, can then be expressed as 



The first factor. 



^a l-(l-7ra)fi<7 



can be interpreted as the transfer efficiency between the curvaton and radi- 
ation. Its maximal value, .^ = 1, is reached when all the energy stored in the 
curvaton is converted into radiation, i.e. when 7^0- = 1, as usually assumed 
in most works on the curvaton. However, if a fraction of the curvaton en- 
ergy goes into species other than radiation, then the transfer efficiency ^ is 
reduced. The second factor, 

(39) 



A-n 



is the familiar coefficient that appears in the literature on the curvaton, which 
coincides with our r only if .^ = 1. 

9 



4.1.2 Second order 

The expressions for the curvature perturbations up to second order are ob- 
tained from the general expression (27-28), using the transfer matrix (32). 
The expression for CDM is relatively simple: 

u = (1 - /c)Cc- + fcC.- + Im - fc) (Cc- - c-f ■ (40) 

The expression for radiation is much more complicated in general, but in the 
limit Xc = 0, which is of interest to us, the radiation perturbation reduces to 



Cr+ — Cr- + - ^a- + — 



2r r 

4r H 



Si- ■ (41) 



In the limit 7j.o- = 1, i-e. ^ = 1, one recovers the usual expression. 

Note that, although fi^ is assumed to be very small, it cannot be neglected 
in the expression for fc [see (25)] because 7co- or Vta^ can be very small, and 
fc can thus take any value between and 1. 

4.2 Initial curvaton perturbation 

We now need to relate the perturbation of the curvaton fluid with the fluctu- 
ations of the curvaton scalar field during inflation. For simplicity, we assume 
here that the potential of the curvaton is quadratic. 

Before its decay, the oscillating curvaton (with mass m ^ H) is described 
by a pressureless, non-interacting fluid with energy density 

Pa = mV' , (42) 

where a is the rms amplitude of the curvaton fleld. Making use of Eq. (9), 
the inhomogeneous energy density of the curvaton can be reexpressed as 

p, = p.e3(C.-5A^) . (43) 

In the post-inflation era where the curvaton is still subdominant, the spa- 
tially flat hypersurfaces are characterized by 5N = (j. (since CDM is also 
subdominant). On such a hypersurface, the curvaton energy density can be 
written as 

p^e^^'^''-'^'-) = p^e"^" =m^{a + 5af . (44) 

Expanding this expression up to second order, and using the conservation 
of 5a /a in a quadratic potential, we obtain 

S, = 2^-f^)\ (45) 

10 



where the initial curvaton field perturbation, (5(T*, is assumed to be Gaussian, 
as would be expected for a weakly coupled field. The curvaton entropy 
perturbation (45) thus contains a linear part S which is Gaussian and a 
second order part which is quadratic in S: 



S^ = S--S\ with S = 2^ 

4 (7* 



(46) 



where the power spectrum of S, generated during infiation, is given by 



{S{k)S{k')) = {2^f -^ Vsik) S{k + k'), Vsik) = -A^ 



(47) 



The subscript * means that the quantity is evaluated at the time when the 
corresponding scale crossed out the Hubble radius during infiation. 

4.3 Primordial adiabatic and isocurvature perturba- 
tions 

For simplicity, we now assume that the post-infiation perturbations for dark 
matter and radiation, i.e. before the curvaton decay, are the same and depend 
only on the infiaton fiuctuations. 



Sc— ^r— Sinf 5 



(48) 



so that there are only two independent degrees of freedom from the infiation- 
ary epoch, ^inf and S. 

Substituting (46) and (48) into (41) and (40) yields 



Sr Sint 2 2g 






S' 



(49) 



and 



S, = ifc-r)S+- 



3/e(l - 2/e 



r(3-8r + — -2- 



^2 



(50) 



In the limit 7^0- = 1, i.e. ^ = 1, one recovers the well-known expression for 

Considering only the linear part of (49), one finds that the power spectrum 
for the primordial adiabatic perturbation ^j. can be expressed as 



Vt 



n., + -w'Ps 



r 
"9 



--1 



:i + A)Pc.„, = 5- -Vs 



(51) 
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where A is defined as the ratio between the curvaton and infiaton contribu- 
tions and S = (1 + A^^)^^ as the ratio between the curvaton contribution and 
the total curvature power spectrum. The hmit A 3> 1, or S ~ 1, corresponds 
to the standard curvaton scenario, where the infiaton perturbation is ig- 
nored. The cases where the infiaton contribution is not negligible correspond 
to the mixed inflaton-curvaton scenario [22]. The curvaton contribution is 
subdominant when A ^ 1, i.e. S ^ 1. 

Let us now turn to the primordial isocurvature perturbation. As can be 
read from the linear part of (50), its power spectrum is given by 

Vs. = {h-rfVs. (52) 

and the correlation between adiabatic and isocurvature fluctuations is 

C ^ -^^^ = Sf S^/^ , Sf ^ sgn(/, - r) . (53) 



In the pure curvaton limit (S ~ 1), adiabatic and isocurvature perturbations 
are either fully correlated, if Sf > 0, or fully anti-correlated, if Ef < 0. In 
the opposite limit (S ^ 1), the correlation vanishes. For intermediate values 
of S, the correlation is only partial, as can be also obtained in multifield 
inflation [35]. 

As combined adiabatic and isocurvature perturbations lead to a distor- 
tion of the acoustic peaks, which depends on their correlation [36], it is in 
principle possible to distinguish, in the observed fluctuations, the adiabatic 
and isocurvature contributions. So far, there is no detection of any isocurva- 
ture component, but only an upper bound on the ratio between isocurvature 
and adiabatic power spectra, which, in our case, is given by 

The observational constraints on a depend on the correlation. Writing a = 
a/(l — a) (note that a ~ a if a is small), the constraints (WMAP-I-BAO-I-SN) 
given in [8] are 

oo < 0.064 (95%CL), a^ < 0.0037 (95%CL) (55) 

respectively for the uncorrelated case and for the fully correlated case ^. 

One sees that the observational constraint a ^ 1 can be satisfied in only 
two cases: 



"'^Our notations differ from those of [8]. Our a corresponds to their a and our fuUy 
correlated Umit corresponds to their fully anti-correlated limit, because their definition of 
the correlation has the opposite sign. 

12 



• \fc — t| <^ r, i.e. a fine-tuning between the two parameters fc and r. 
This includes the case /c = 1 with r ~ 1, considered in [17]. 

• S <^ 1, i.e. the curvaton contribution to the observed power spectrum 
is very smaU. 

4.4 Non-Gaussianities 

Let us now examine the amphtude of the non-Gaussianities that can be 
generated in our model. We recall that our observable quantities ( and S are 
of the form 

C = Chif + ziS + -Z2S^, S = siS + -S2S^, (56) 

where Cinf and S are two independent Gaussian fields and where the coeffi- 
cients can be read explicitly from (49) and (50). 

Applying the general results of the Appendix to the present situation, 
we can easily compute for our model the "reduced" angular bispectrum, 
which is of direct interest for a comparison with CMB observations and which 
generalizes the analysis of [37] in the purely adiabatic case. Specializing (112) 
to our case, one finds 

I,J,K 

with 



W.is = 3 }^ 6^f / r2rfr/3fjr)A^(r)/35(r) (57) 



&^f = Arf2)Ar({)iV(f), (58) 

where A/^L = ^2, N^^. = S2, A^^-, = zi, A'^n = si, respectively, and 

/3/(r) ^^J edkMkr)gl{k), ^l{r) ^^J edkMkr)g^{k)Ps{k) , 

(59) 
where the gl{k) denote the adiabatic (J = () and isocurvature (/ = S) 
transfer functions. Because of the symmetry under exchange of the last two 
indices, the reduced bispectrum contains six different contributions, whose 
amplitude is parametrized by the six coefficients b^j^ . 

In order to compare these coefficients with the usual parameter /jvl de- 
fined in the purely adiabatic case, one must recall that /nl is proportional to 
the bispectrum of ( divided by the square of its power spectrum. By noting 
that the /3/(r) introduced in (59) involve Pg, this implies that the analogs 
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of Jml can be defined by dividing the coefficient 6^^. by the square of the 
ratio Pc^l Pg = zfS"'^. We thus introduce the parameters 

fI,JK _ pI,JK _ ^ hI,JK ff^r,\ 

Jnl — r^JNL — 4 "nl ■> y^^i 
explicitly given by the expressions 

f CCC _ 22 ■='2 f CCS _ siZ2 'z:2 K,SS _ sfz2 ^2 /ei \ 

JNL ~ I? " ' JNL ~ 2? " ' JNL ~ ~7r ^ ' '^ '' 



fS,QC, _ S2 -^2 fS,i;S _ SIS2 -^2 fS,SS _ sfs2 ^2 (cn\ 

JNL ~ i? " ' JNL — ~TS- ^ ' JNL ~ ^F " ' v" '' 



In the absence of isocurvature perturbations, the above non-linear pa- 
rameters vanish except the first one, yielding 



5 


/ 3 2 




r 




\ 


-4- 
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\2r e 




e 



/^r = ^ £: + 7-4-7^ H% (63) 



which exactly coincides with the familiar parameter /atl. The amplitude 
of the non-gaussianities is determined by the three parameters r, ^ and S 
(note that one recovers the usual prediction of the pure curvaton scenario 
for ^ = 1 and S = 1), which take values between and 1. The present con- 
straints on J'nl, calculated from WMAP data by assuming purely adiabatic 
perturbations, are [8]: 

-10 < /Sr'^ < 74. (64) 

A sufficiently small r, or even ^, leads to a significant non-Gaussianity from 
the adiabatic component, whereas a small H tends to suppress it. 

If isocurvature modes are present, however, the five other terms in the 
reduced bispectrum (57) will also contribute in general. Interestingly, the 
six functions on the right hand side of (57) have distinct dependences on the 
li, because they involve different combinations of the adiabatic and isocur- 
vature transfer functions. Therefore, this allows in principle to measure, or 
constrain, independently the corresponding six non-linear parameters from 
the CMB data. The precise determination of constraints on the f^i^ is 
beyond the scope of the present work, but since all the functions multiply- 
ing the bj^i^ in (57) are of similar amplitude, one can a priori expect the 

rI,JI' 
'NL 



constraints on the f^j^ to be of the same order of magnitude as those on 



Jnl ■ 



^Observational constraints on isocurvature non-Gaussianities are given in [19] , for an 
isocurvature perturbation of the form S = Sl + f^L ^l ' where Sl is Gaussian. Their 
non-hnear parameter /j^^° is related to ours according to /^'^ — 2fj^'^^'a'^, fj^'f^ = 
2/^^°^a3/2|^| j^j^j fSXC ^ 2/^^°)aC2, where C is the correlation defined in (53). 
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Let us now explore the amplitude of the non-linear parameters in our 
model. First of all, let us stress that finding significant non-Gaussianities 
(typically /jvl ~ 10 — 100) requires, in all cases, a small denominator zi, i.e. 
r <^ 1, which will thus be assumed below. Second, it is worth noting that all 
the coefficients are related via the two rules 

Inl = —Inl (^1). Inl = —Inl > (^2) . (65) 

Therefore, the hierarchy between the parameters can be deduced from the 
value of the first order ratio 

— = 3 ( — -1 I =£/W^, £/ = sgn(/e-r), (66) 

where we have used (54), as well as the second order ratio S2/-225 which is a 
more complicated expression in general. 

We now consider successively the two limits for which the isocurvature 
bound is satisfied. 

4.4.1 Limit |/c — ^| ^ r, with S ~ 1 (pure curvaton scenario) 

In this case, the isocurvature-adiabatic ratio a must satisfy the observational 
constraint a ~ Oi < 0.0037, since we are in the fully correlated case. The 
relevant ratios are given here by 

-f^-^'^-f) (67) 

l + 2f(2-f)/3 ^ ' 

where we have taken the limit r = ,^ f ^ 1 (although r cannot be smaller than 
10~^, to be compatible with observational constraint on Jml)- If ^ is small 
because f <^ 1, then the denominator in the expression for S2/ Z2 reduces to 
1. However, if ^ <^ 1 while f is of order 1, the full expression for S2/ Z2 is 
needed. 

The value of the first ratio implies that, with respect to /^'^ , the coeffi- 
cients /^'2^ and fj^j^ are suppressed with factors y/a and a, respectively, ac- 
cording to (Rl). Analogously the coefficients /^'^ and f^j^ are suppressed, 
respectively with factors ^/a and a, with respect to /Vl • By contrast, using 
(R2), one sees that fj^jj could be of the same order of magnitude as /^f, , if 
f ~ 1, or suppressed if f is small. 

To conclude, in the pure curvaton scenario, it is possible to satisfy the 
isocurvature constraint and to get measurable non-gaussianities only by as- 
suming a fine-tuning between /c and r at the percent level. In this situation, 
only the purely adiabatic parameter is significant, while the other parameters 
are suppressed, with increasing powers of a. 
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4.4.2 Limit S < 1 

In this limit a must satisfy tlie constraint a ~ oq < 0.064 (uncorrelated 
case). 

In the regime /c ^ r ^ 1, one finds that both ratios s^j Z\ and s^j z^ re- 
duce to (-3), independently of the value off. Therefore, the relation between 
the non-linear parameters is simply 

/^f=^^, /iL"" ^ (-3/^/Jf (/c«r«l) (68) 

where Is is the number of S among the indices (J, JK\ This is the result 
obtained in [17] for fc = 0. For a close to its present upper bound, one sees 
that detectable non-Gaussianity can be generated with r ~ 10~^. 

By contrast, in the regime fc ^ r, the purely adiabatic coefficient is 
strongly suppressed since 

However, the other coefficients are now enhanced with respect to the purely 
adiabatic coefficient, via the large factors 

-~3^, £i^3^(l-2/,). (70) 

Zi r Z2 r 

where, for simplicity, we have assumed f ^ 1 (the other possibility .^ ^ 1 
yields a more complicated expression for the second ratio, with a dependence 
on f). The dominant term is therefore the purely isocurvature term 

S'-^^m^ (71) 

If /c ~ 1, this purely isocurvature non-Gaussianity, although enhanced with 
respect to all the other contributions, remains negligible since it is suppressed 
by the very small factor a^. This was the conclusion reached in [17] (for 

fc = 1). 

However, we now see that this suppression can be compensated if fc is 
smaller than a"^. The purely isocurvature parameter and the other ones are 
then given by 

/^f^^, fYL^" ^ (Jj-y^ fjf (r«/c«l) (72) 

where /^ is the number of ( among the three indices. One can notice that the 
amplitude of the purely isocurvature non-Gaussianity does not depend on the 
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parameter r, but only on a and /c. For instance, with a = 0.05 which satisfies 



the current observational bound, a value /c = 10 yields /^ 



s,ss 

NL 



100. In 



such a scenario, one gets observable non-Gaussianity that conies essentially 
from isocurvature modes, even if the latter are subdominant in the power 
spectrum. 

5 Scenario with two curvatons 

We now apply our formalism to the models where two curvatons are present 
in the early Universe (see e.g. [25, 26, 27]). The curvaton a will be assumed 
to decay first, followed later by the curvaton denoted x- 



5.1 First order 

At linear order, the decay of the first curvaton can be characterized by the 
transfer matrix 



Tr 



[1] 



i-fci u 

\ / 



(73) 



where the order of the species is (r, c, XjO"), while the decay of the second 
curvaton is characterized by the transfer matrix 



Tr 



[2] 



/I 



V 



-a;^2 


Xc2 


Xr2 — Xc2 


M 





l-/c2 


fc2 


























0/ 



(74) 



In the above matrices, the definitions of the parameters are analogous to 
the definitions introduced in (32), i.e. x^i = fri/Cli, Xd = QciXri/4:, x^i = 
^x^-'^ri/^:, etc, and the indices 1 and 2 refer respectively to the first and 
second decays. We have also allowed the possibility that the first curvaton a 
decays into the second curvaton x, hence the presence of the parameter /^i. 
The expression of the perturbations for radiation and CDM, after the 
two transitions, are expressed in terms of the initial perturbations (bo via 
the product of the two transfer matrices given above, i.e. 



Ca = 5^(T[2] -^[1])^ CbO- 



(75) 
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At first order, the radiation curvature perturbation, after tlie second cur- 
vaton decay, reads 



Cr — CrO + Z(J 5*0-0 + % 5*^0 + ^c 'S'cO, (76) 



witli 



"iZa = {l-Xr2){Xrl-Xcl-X^i) + fclXc2 + fxl{Xr2-Xc2), (77) 

3z^ = (1 - fxi){xr2 - Xc2) + (1 - a;r2)a;xi, (78) 

3Zc = (1 - fcl)Xc2 + (1 - Xr2)Xcl ■ (79) 

Combining tliis expression witli tliat of tlie CDM curvature perturbation, 
according to (11), we find tliat tlie CDM entropy perturbation is given by 

Sc = So- 5*00 + Sx ^xo + ^c Sco, (80) 

witli 

So = -3^o + /cl(l-/c2) + /c2/xl, (81) 

s^ = -3z^ + U{l-fxi), (82) 

s, = -3^, + (l-/,i)(l-/c2). (83) 

For simplicity, we will restrict ourselves, from now on, to the case where 

Sco = 0. 

Defining A as the ratio between the two curvaton power spectra, such 
that 

Ps^o ^ AP5.0, (84) 

one easily finds that the ratio between the isocurvature and the adiabatic 
spectra is given by 

where A^ and Ao are defined as in (51), i.e. 

n = ^H, + 4^5.0 + 4^5X0 = (1 + Ao + Xx)'PCrO- (86) 

The correlation between Cr and Sc can be expressed as 

[4 + A4)(4 + A4) 
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The observational constraints on a impose that at least one of the follow- 
ing conditions must be satisfied: 

S < 1 or s^ + AsJ < 2^ + kz\ . (88) 

The first possibility, S <^ 1, corresponds to a power spectrum dominated by 
the inflaton, whereas the second possibility requires special cancellations in 
(81-82) so that s^- and s^ are suppressed. 

5.2 Second order 

We now consider the perturbations up to the second order, in order to com- 
pute the non-Gaussianities. First, let us decompose the curvaton entropy 
perturbations as in (46), so that 

Scto = So- — jS^ S^o = S^ — jS^, (89) 

where S„ and S^ are two independent Gaussian quantities. 

The radiation curvature perturbation and the dark matter entropy per- 
turbation after the second decay, up to second order, are given in our notation 

by 

^ ^ 1 - 1 - 

Sc = ScrScT + s^S^ + Sf^^StjS^ + -SfjcrS^ + -Sxx^^ (91) 

where the coefficients z„, z^, s„ and s^ have already been defined in (77- 
78) and (81-82), respectively. We do not give explicitly the full expressions 
for the second order coefficients because they are very lengthy, but they are 
straightforward to compute by using the general expressions (27-28) with the 
transfer matrices (73-74). 

Let us calculate the reduced bispectrum by using the general expression 
given in the Appendix. In our model, ignoring the inflaton which does not 
produce signiflcant non-Gaussianities, the relevant power spectra are inde- 
pendent so that 

P'\k) = ( J ° ) P,^ , (92) 

where we furthermore assume that A is strictly independent of k (this is 
indeed the case if the masses of both curvatons are negligible with respect to 
H during inflation). 
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As a consequence, the reduced bispectrum can be reduced to the same 
expression as that already given in Eq. (57) with 



/3/(r)^- I kHkji{kr)gl{k)P^^{k) 



(93) 



and the six parameters 



,I,JK 
-'NL 



Ni^NiN^ + AiV^ {N^N^ + N'^N^) + K^N'^^N'^N^ , (94) 

where the coefficients A^^^, which are defined as in (105), can be read off 
directly from (90) and (91). In complete analogy with the model with one 
curvaton, to be compared with the standard J'nl, these coefficients must be 
divided by the square of the ratio Pf^/Ps^^ = {z^ + Az^)/E, hence: 

2 



71, J K 

Jnl 



zl + Kzl 



,I,JK 
-'NL ■ 



(95) 



The six non-linearity coefficients are thus given by 



fCCC 



JNL — 


u 


+ kzl 


Jnl — 


/ 


'z: 


V4 


+ A4 


K,SS _ 


/ 


'Z; 


Jnl 


VI 


+ Azl 


Jnl — 


/ 


'z 


V4 


+ Azl 


Jnl — 


/ 


^ 


1 4 


+ A4 


fS,SS 


/ 


^ 


Jnl 


\4 


+ Kzl 



^aa^iy ~T ^i\-Zrj-^ZfjZy^ + i\ ^^xX^xJ ' 

^aaZa^a ~r -li- Z(JX\Za^x ~r ^x '^1 ~^ XX X xj ' 

^aa^a + li\ Zfy-^SfjS-^ + A ^^X-^xJ ' 

°aaZ^ \ ^■^'■^CTX cr^x ' XX xJ ' 



Saa^aSa ~r ^^^axy^o-^x ' "^X^""/ ' '^XX'^X^xJ ' 



s^^s^ + 2As^^s^s^ + A^s^^sJ] . (96) 



In the following we analyze explicitly some limiting cases. 

5.3 Various limits 

We now explore the parameter space, in order to see whether it is possible 
to obtain significant non-Gaussianities. 

Let us first mention that we have checked that our results agree with 
those of [27] in the limit where the curvatons decay only into radiation (i.e. 
/ci = /c2 = /xi = 0); the dark matter abundance is neglected (i.e. Xd = 
Xc2 = 0) and the inflaton contribution is ignored (i.e. S = 1). 
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5.3.1 Limit A < 1 

In this limit wliere tlie contributions from the second curvaton are neghgible, 
one finds 

«--4 /?i-S^^, (97) 

while the other five non-hnear coefficients can be deduced from fj^l according 
to the relations 

f ^-^5 _ S^ rI,JC rS,IJ _ foa K,IJ /qs^ 

JNL — JNL JNL — JNL ■ \^°) 

The quantity a is constrained by observations to be small, which requires 
either S ^ 1 or |so.| ^ \z„\. 

First possibility: S ^ 1, while \za\ ~ \sa\- 

This leads to a suppression of all the non-Gaussianity coefficients by a 
factor 'E?. However, the coefficients /^'^ can still be significant if the ratio 
z„^l z\ can compensate the 'E? suppression (similarly for the f^^ if s„^/ z^. 
compensates the E? suppression). 

Let us consider a specific example, with the simplifying assumptions 

Xci = Xc2 = fxi = a;^! = , (99) 

that is, we neglect the energy fraction of dark matter and assume that the 
curvaton a does not decay into x ^-iid that x is subdominant when a decays. 
Under these assumptions, z^ = Xri{l—Xr2)/3 and we further assume /d <^ z^ 
so that So- — —3Z(j. In the two hmits Xri = riC,i ^ 1 and (1 — Xr2) ^ 1, Za 
is small and the adiabatic non-Gaussianity behaves as 



fCCC _ L 

JNL ~ 



Xr2 



K,CC I ^r2 (3 pr^_^^ 



(100) 



where /^ii2 correspond to single-curvaton coefficient, equation (63), but 
calculated with the parameters ^1,2 and Xri,r2 respectively. 

If we assume Xr2 ^ 1, the above equation corresponds to the single- 
curvaton result (68). The other coefficients also follow the relations given in 
(68), since So-a/zaa = —3 with the assumptions (99) and /d <^ 1, and are 
thus of comparable magnitude. 
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Second possibility: \scr\ <^ \za\ 

When a small a is the consequence of \s„\ <^ \z„\, one sees from the first 
relation in (98) that all the fj^^ are strongly suppressed with respect to Z^^. • 
However, the two coefficients fj^}^ can still be important if \zcr\ is sufficiently 
small. By examining (77) and (81), one sees that getting \s^\ <C |^o-| ^ 1 
requires some fine-tuning between the coefficients, which we now discuss. 

In order to get \zcr\ ^ 1, the first possibility is that the first curvaton is 
subdominant, i.e. x^i = 0{e), where e is some small number (we neglect Xc2 
which must be small because we are deep in the radiation era), which then 
requires either Xr2 = 0{e) or /^i = C(e). The second possibility is that the 
second curvaton dominates at decay, i.e. Xr2 = 1 — ^^(e)? which also requires 
that /^i = 0{e). Then, to obtain |so-| ^ \za\, the terms of the right hand side 
of (81), which are of order e must compensate each other so that their sum 
is at most of order 0{ae), which necessitates some special relation between 
the Ja and the xa- 

If we consider again the assumptions (99) and neglect fc2, one finds that 
the fine-tuning condition to get |scr| ^ \za\ is 

fci-xriil-xr2)<0{ae). (101) 

The adiabatic parameter /^f, is given in equation (100), with now S ~ 1, 
and its value is of order 10 when e ~ Xri{l — Xr2) ~ 0.1. Since Saa/zaa — 
-3 + 0{fci/xri{l - Xr2)), we also have f^i^ ~ f^'l^. 

Note that a significant non-Gaussianity generated by a dominant curvaton 
{Xr2 = 1 — C^(e)) has already been pointed out in [27], but we see here 
that satisfying the isocurvature bound requires additional constraints on the 
branching ratios of the curvatons. 

5.3.2 Limit A > 1 

In this limit, one obtains 

^ ^ fx KCC _ ^^2 £xx fi,JS _ £x ri,JC rS,iJ _ £xx K,iJ 

" "" ^2 ' JNL " 2 ' JNL — JNL ' J NL — J NL ' 

^X ^X X ^XX 

(102) 
By comparing with (97) and (98), one sees that the analysis is analogous to 
the previous case, by replacing z^-, z^^^, s„ and s^^ by z^, z^^, s^ and s^^, 
respectively. 

When the curvaton contribution to the power spectrum is not negligi- 
ble, significant non-Gaussianity, while satisfying the isocurvature bound, is 
obtained when \s^\ <C \z^\ <^ 1. This constraint is satisfied if one assumes 
f^i = 1 — C(e), which means that the second curvaton is created mainly by 
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the decay of the first, while Xr2 = 1 — 0{e), x^i < (9(e) and /c2 = 1 — C^(e)- 
Other possibilities exist but require some fine-tuning between the parameters, 
in analogy with the previous analysis in the case A ^ 1. 

5.3.3 Intermediate values of A 

In this case, one must satisfy simultaneously the constraints |so-| ^ ko-l and 
kxl ^ kxl' ^^^ ^° ^^^ isocurvature bound. The relative strength of the 
different j'nl coefficients cannot be expressed in such a simple form as in 
(98), but it will be determined again by the ratios Sa/z^, s^/z^, Saaj^aa and 

In order to get also a significant non-Gaussianity, we look for parameter 
values such that 

z,,z^^O{^, s,,s^<0{at). (103) 

These constraints can be satisfied by fine-tuning the parameters. Solving 
So- ~ and s^ — '-' ^^^ ^^^ ^^° parameters /d and fc2 yields 

^^^ ^ (x,.-x..)(l-/xi)-:..i _ ^^^ ^ ^^^ - -.. + 1^ -XI . (104) 

J- - Jxi - 2:xi i - /^i 

The observational constraint on the isocurvature power spectrum is satisfied 
if these two fine-tuning relations hold simultaneously, at the level Oiae). 
Using these relations, one finds interesting non-Gaussianity for the following 
set of parameters: x^x = 0{e), Xr2 = 0{e), x^i = 0{ae), /d = Xri — Xd + 
0{ae), fc2 = Xr2 + 0{at), with negligible values for Xc2- In this scenario, 
both curvatons are subdominant at their decay and the fraction of produced 
dark matter is fine-tuned. 

6 Conclusions 

In this work, we have introduced a systematic treatment in order to com- 
pute the evolution of linear and non-linear cosmological perturbations in a 
cosmological transition due to the decay of some particle species. Our main 
results can be summarized as follows. 

At the linear level, the evolution of all individual curvature perturbations 
can be expressed in terms of a transfer matrix, whose coefficients depend on 
background quantities, such as the relative abundances of the fluids at the 
decay, their equation of state parameters and the relative decay branching 
ratios [see Eqs (21-25)]. At the non-linear level, the post-decay curvature 
perturbations can also be given in terms of the pre-decay perturbations quite 
generally, and we have presented explicitly these relations at second order [see 
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Eqs (27-28)]. We have then apphed our general formahsm to two specific 
examples. 

The first example is the mixed curvaton-infiaton scenario in which we 
allow the dark matter to be created both before and during the curvaton 
decay. We find, in particular, the remarkable result that it is possible to 
obtain isocurvature dominated non-Gaussianities with, as required by the 
CMB measurements, an adiabatic dominated power spectrum. 

In the second example, we have studied scenarios with several curvaton- 
like fields and obtained results that generalize previous works on two-curvaton 
scenarios by taking into account the various decay products of the curvatons. 
We have explored the parameter space to see whether it is possible to find 
significant non-Gaussianity while satisfying the isocurvature bound in the 
power spectrum. We have found that several such regions exist, but often at 
the price of a fine-tuning between the parameters. 

In the presence of isocurvature modes, which can be correlated with the 
adiabatic modes, non-Gaussianity of the local type is much richer than in 
the purely adiabatic case and we have shown that the angular bispectrum is 
the sum of six different contributions. As a consequence, in addition to the 
traditional /tvl (adiabatic) parameter, we have identified five new non-linear 
parameters that must be taken into account: one purely isocurvature param- 
eter and four correlated parameters. We have computed these parameters in 
the two models we have investigated. 

Beyond the two examples considered in this work, our formalism can 
be used as a general toolbox to study systematically the cosmological con- 
straints, arising from linear perturbations and from non-Gaussianities, for 
particle physics models and their associated cosmological scenarios. 
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A Angular bispectrum 

We consider several observable quantities X^ , like ( and Sc, which depend 
on "primordial" scalar fields 0", whose perturbations are generated during 
infiation. Up to second order, one can formally write 

X^ = Ar„V + iAr>V + ... (105) 
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We assume that the (p"" are Gaussian random fields, with the two-point cor- 
relation functions 

{<P''{k)<p\k')) = {27T'f P^\k) 5(k + k') . (106) 

We then define the bispectra of the X^ by 

{Xlxixg) = {2nfS{E,h)B''''{k,, k,, k^) . (107) 

Substituting the decomposition (105) into the left hand side, and using (106), 
one finds 

i?"^(A;i,A;2,A;3) = NiNiN^^P'^%k,)P'>%k2) + NiXiN^ P'^\k,)P'\k^) 

+Xi,XiX^P^\k2)P'\k,). (108) 

As shown in [37], the angular bispectrum can be expressed in terms of 
the "reduced bispectrum" bi-^^i^i^, according to the expression 

where 

QT.l.rr' ^ jd'nYi.^Ai^) Yi,^,{n) Yi^^^{n) (110) 

is the Gaunt integral. 

The next step is to express the reduced bispectrum in terms of the gener- 
alized bispectra B^'^^ , using the fact that the observable quantity is related 
to the initial perturbations X^ via a transfer function g{{k), so that 

^^m = M-^y J -0^, fE^'(^)^'(^)) ^™(^)- (111) 

Substituting the above expression in the left hand side of (109), one finally 
obtains 

6,,y3 = 3 J2 NLNiNf / r^rfr^fjr)Af ^(r)/35^^(r), (112) 

I,J,K -^0 

with 

/3/(r) = ^ y k'dkjiikr)glik), Pl'^^^r) ^^J k^dkji{kr)gl{k)P-\k) . 

(113) 
Note that the "reduced" bispectrum is symmetric with respect to permuta- 
tions of the indices /i, I2 and l^ (we use the standard notation: {hhh) = 
[hkh + 5 perms] /3!). 
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In the simplest case, one considers only the adiabatic mode, ( (or the 
gravitational potential $), which is assumed to depend on a single "primor- 
dial" Gaussian field. In this case, where both the indices / and a take a 
single value, one recovers immediately the familiar result of [37]. Our gen- 
eral expression also includes the particular situation considered in [19], where 
C = + (3/5)/ArL0^ and 3 = 1] + /jvl^^, and tj being Gaussian variables. 
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